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Abstract 

We reexamine the problem of having nonconservative equations of mo- 
tion arise from the use of a variational principle. In particular, a formalism 
is developed that allows the inclusion of fractional derivatives. This is done 
within the Lagrangian framework by treating the action as a Volterra series. 
It is then possible to derive two equations of motion, one of these is an 
advanced equation and the other is retarded. 



1 Introduction 

The problem of having a dissipation term q arise in the equations of motion for 
a system has a long history. Bauer J21 showed that "the equations of motion of 
a dissipative linear dynamical system with constant coefficients are not given by 
a variational principle". There are loopholes in Bauer's proof, however. One of 
these is to allow for additional equations of motion to arise. This method was 
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employed by Bateman [ 1 1. He used the Lagrangian 

C 

L = raxy + — (xy — xy) (1) 

which gives the equations of motion 

rax + Cx = my — Cy = (2) 

Bateman's method is not very general, so we look for other methods to model 
nonconservative systems. 

Caldeira and Leggett [6| suggest recognizing that a dissipative system is cou- 
pled to an environment. The environment is modelled as a collection of harmonic 
oscillators which results in the Lagrangian 

l = ^e-v( q ) + ±^e n -^( qn - q r} (3) 

where q is the system's coordinate and the g n 's are the environment's coordinates. 
While the system by itself is nonconservative, the system plus environment is con- 
servative. This procedure does allow the introduction of very general dissipation 
terms into the system's equation of motion. However, the microscopic modelling 
of the environment makes © much more complex than, say, O. 

In order to overcome the difficulties of the above two procedures, Riewe 
examined using fractional derivatives in the Lagrangians lfT2l 031. This method 
takes advantage of another loophole in Bauer's proof. Namely, Bauer assumed 
that all derivatives were integer ordered. Riewe's method has the advantage of 
not introducing extra coordinates as in and ©. However, it ultimately results 
in noncausal equations of motion. A rather ad hoc procedure of replacing anti- 
causal with causal operators needs to be used at the end in order to arrive at causal 
equations of motion. We will present a method that can be used within Riewe's 
formalism that avoids this situation. 

We propose here a new method of using a variational principle to derive 
nonconservative equations of motion. Our method is closely related to Riewe's in 
that we use fractional operators. However, we treat these operators as kernels in a 
Volterra series. We show that Riewe's formalism can be derived by using certain 
types of symmetric kernels in the series expansion. A simple modification of the 
kernels will result in two equations of motion for a system. One of these equations 
is advanced while the other is retarded, similar to ©• 
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Our paper is organized as follows. In Section|2]we review fractional integration 
and differentiation. Riewe's formalism is briefly examined in Section|3] We then 
give a brief overview of Volterra series in SectionHJbefore examining our fractional 
Lagrangian mechanics in Section Section ]§\ examines the nonconservative 
harmonic oscillator in a different way than the traditional variational methods. A 
discussion of some related concepts and future research follows in Section|7] 



2 Fractional Integration and Differentiation 

Fractional integrals and derivatives are generalizations of their usual integer or- 
dered operations. To start developing the theory, let us first write down Cauchy's 
integral formula 

™ = rR />><<- T >"-'* (4) 

where n > is an integer, T(n) is the gamma function, and a < t. Equation © is 
a convolution of f(t) and the function 

$+U) . = J r n) 1 1 > u (5) 

if we set f(t) =0 for t < a. So we can rewrite © as 

.!?[/] = /(*)*<#(*) (6) 
where * is the convolution operation defined by 

/oo 
g(r)h(t-r)dr (7) 

Equation © will be our stepping stone to generalizing the integer ordered opera- 
tions to fractional order. 

The above procedure works so well for the integers n > 0, we want to consider 
extending it to any real a > 0. This is obviously possible, so we let 

a I t *[/] = /(*)*<&+(*) (8) 
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be the left fractional integral [LFI] of f(t) of order a > 0. Everything works fine 
until we consider the case a = 0. We reasonably expect that 



aVlf] = /(*) (9) 

but, it is not immediately obvious that the integral in © is not divergent. Also, for 
— 1 < a < 0, the integral is obviously divergent. It is apparent that treating f(t) 
and $+ (t) as regular functions will not be sufficiently general for our purposes. 
Instead we will consider them to be distributions, or generalized functions. [We 
note that there are other ways to generalize integer ordered derivatives to frac- 
tional order iTTOl . We will work solely with the generalized function approach to 
interpolate between the integer ordered integrals and derivatives.] 

The first order of business is to define the convolution operation for distribu- 
tions. Let k(t) = g(t) * h{t) and cp(t) be a test function. Then 



(k,<p) := / k(t)<p(t)dt 

J g(Oh(t-Z)dt}<p(t)dt (10) 



Equation (fTOb is meaningful as long as either g(t) or h(t) has bounded support or, 
g(t) and h(t) are bounded on the same side [e.g., g(t) = for t <ti and h(t) = 
for t < t 2 ]. We will always assume that one of these situations is the case. From 
(fTOl. it can be seen that the generalization of © is 

(g*h,<p) = {g(t),(h(T)Mt + T))) (ID 
The convolution operation has the properties 

g * h = h* g (12) 
f*(g*h) = (f*g)*h (13) 
D(g*h) = (Dg) * h = g * (Dh) (14) 

where D(-) is the generalized derivative. Remember that the relationship between 
the generalized and classical derivatives, beginning at t — a, is given by ifTUl 

n— 1 

D n f = + [D n ~ fc-1 tf(t - a)] f ik \a) (15) 

fc=0 
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where is the classical derivative. 

Considering as a generalized function allows us to extend © to any a, 

where the convolution operation is defined as in dTTb . For a < 0, this will define 
the left fractional derivative [LFD] as 

a D-«m := a I?[f] 

= f(t)*®t(t) (i6) 

In the sequel, we will find it easier to assume a > and use the notation 

8 D?[/] = /(*)**+,(*) Dr a [/] = J?\f\ (17) 

Also, for reasons that will become apparent shortly, we will often set f(t) = for 
t < a and t > b, where a < b. We do not want any resulting discontinuities in 
f(t)&tt = b to affect the LFDs. So t must be restricted to the interval a < t < b 
in the LFDs. It would perhaps be better to write (fTTT) as 

a -D?-[f] = 7^ f f{r){t-r)-^dr (18) 



r(-a) 

To avoid cluttering our notation, we will continue to use the notation in (fT7l with 
the understanding that it formally means (fTSl . 

The distributions <E>+ (i) have been well studied J7J [TU||. Their two most 
important properties are 

$+(*) = £T"5(t + ) (19) 
for any integer n, and, for any (3 and 7, 

a) = $J +7 (t-a) (20) 

Equation d2"0l implies 

a Df[ a D7[/]] = a Df + ^[/] (21) 
aD^ml = / (22) 



Now let < n - 1 < a < n. Then, using <H2) - ([HI) and <Q]3) and (EOT) , we 

have 

aPt if] = f(t) 

= (D n f(t))*^- a (t) (23) 
= D n (f(t) * <&+_ a (t)) (24) 
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Equations d2"3t and d24t> are the distributional forms of the Caputo and Riemann- 
Liouville fractional derivative, respectively |10|. In the standard definitions of 
these derivatives, D n is replaced with (d/dt) n . 

In addition to the left fractional operations, we can also define right fractional 
operations. If we set /(£) = for t > b and define 

the right fractional operations are defined by 

t D?[/] := (26) 

Most of the above observations for the left fractional operations also hold for the 
right ones. However, (fT9"l) needs to be replaced with 

$;(£) = (-l) n D- n 5(t-) (27) 

for any integer n. When /(£) = for t < a and £ > 6, we do not allow any 
resulting discontinuities in /(£) at £ = a to affect the RFDs. Similar to the case 
for the LFDs, we will take (|2"6T) as meaning 

= TT7~ \ f(r)(r-t)-^ (28) 
1 (—a) ./t+ 

though we will continue to use the notation in (l2"6l) . 

Note that for the left operations, the "left" integration limit a determines the 
allowable functions in the operation a D" [/]. Namely, /(£) must vanish for t < a. 
Also, a Df [f] is a function of a and £ and, a functional of /(£). Similar comments 
hold for the right operations. Here, the "right" integration limit b means /(£) = 
for £ > b. Now let /(£) be compactly supported on the interval [a,b]. Then 
a D™ [/] = whenever £ < a. However, a Df [/] does not generally vanish for 
£ > a. Thus, the left operations are causal or retarded. Conversely, t ~Df; [/] = 
whenever £ > b but, generally, t D% [f] ^ for £ < b. Hence, the right operations 
are anti-causal or advanced. 

Our fractional derivatives satisfy an integration by parts formula. First, assume 
that /(£) = for £ < a and g(t) =0 for £ > b. Then, for any (5, 

(g(T)*i(r-t)f{t),<p(t,T)) = (js(r)*j(t-T)f{t)Mt,T)) (29) 
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Hence, 



{g(*t*f),<p) = ((g**p)f,<p) (30) 



or 



J a T>?[f]g(t)dt = j t T>l\g]f{t)dt (31) 

We note that Riewe's derivation of an integration by parts formula |[T3l Equa- 
tion (16)] is flawed on two points. First, the boundary conditions are generally 
fractional, not integer, ordered. Also, Riewe incorrectly exchanges the classical 
Caputo derivative r d2"31 with D n replaced with (d/dt) n ] for the Riemann-Liouville 
derivative in (|2"?1) . Fortunately, when vanishing boundary conditions are assumed, 
these defects are inconsequential. Also notice that (|3"T1) implies that any integration 
by parts inherently introduces time reversal. 

When we examine Riewe's fractional mechanics in Section |3J (|3"T1) will lead to 
equations of the form 

^ *($+*/) = t mipt\f\] 

%p *$+)*/ (32) 



The difficulty with d3*2l is that neither are $t(t) or compactly supported, 

generally, nor are they bounded on the same side. So we need to make sense of 
the convolution in d3*2l . To give meaning to the convolution, let us note that the 
Fourier transform of $^ (t) is given by Q 



exp[sgn(u;)i/?7r/2] 



and for §g (t) 



' lour 



[Note that (1331 and d3~3b imply that, up to a sign, the fractional derivatives go to 
the integer ordered derivatives when (3 is an integer.] Then, 

$-(*)* $+(t) ^ (35) 
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Now, 



2cos(/37r)r(2 / 3) 



T I 1-2/3 
< > \Uj\ 1 



(36) 



We will define 



2cos(/57r)r(2/?) 



and let 



*j(t)**+(t)*f(t) s * v (t)*f(t) 



(37) 



(38) 



for any /3 where /(i) = for t < a and t > b. We call (l38l) a Feller fractional 
derivative [FFD] [ 1 1 1 and write this as 



(39) 



Note that, for n an integer, 



for < t < T, but 



tn n if] = {-i) n f {2n \t) 



(40) 



(41) 



Some care is needed when using the FFDs. Formally we have set f(t) = for 
t > a and t < b. However, the LFD only acts on the resulting discontinuities that 
may be present in f{t) at t = a, not at t = b. Conversely, the RFD acts on the 
discontinuities at t = b, not t = a. It is perhaps better to write (l37l) as 



# 2j s(*) 

Then (l39l) can be written as 



1 



2 cos(/?tt) 



2/3 ^ 



(42) 



2 cos(/?vr) 



{ Df [/]} 



(43) 
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We will take (09"1) as implying d4*5t . 

In general, the fractional derivatives are nonlocal in time. That is, they have 
a "memory". For integer ordered LFDs and RFDs, this memory disappears [i.e., 
they are "amnesiac"] and they act locally in time. Even integer ordered FFDs are 
also amnesiac since the kernels ^ r _ 2n (i) equal, up to a sign, $i 2n (t) and $I 2 n(^) 
in this case. All of the fractional derivatives have a fading memory, however 0. 
That is, they are affected more by the recent past and/or future than the distant past 
and/or future. 



3 Riewe's Fractional Lagrangian Mechanics 



Here we examine Riewe's fractional mechanics lfT^l[T31 . restricting our attention 
to Lagrangian mechanics with Lagrangians of the form 



y U) 2 + f u a ) 2 -n?) 



L{q, a q?, a ql) = - („<//)" - ~ L<l? T - V(<i) (44) 
where q is our [generalized] coordinate, C is a constant, < a < 1 and, 

alt ■= aV"[q] (45) 

We define the action associated with (EH) by 

S[q] := t Ldt (46) 

J a 

Let us consider perturbations rj(t) of q(t) where r](t) vanishes for t < a and t > b 
but is otherwise arbitrary. Then, 

S3[q] = 5 f Ldt 

J a 

~L(q + V , a q? + Jfi, a q\ + A) ~ HQ, a Qt, ati)} dt (47) 
Expanding the perturbed Lagrangian in (13Tb 

L(q + V, a q? + Jit, a4 + aVl) = Hq, X, + 

9L dL a dL , ia 
^ V+ 9^ X+ 0^^ <48) 



and using (B8b in ( R7b . we have 



77 + 



<9L 



<9g d a qf 



f b (dL ^ 



dL 



dL 



Ail. 



dt 



(49) 



where we used (|3"TI in going to the second equality. 



Hamilton's principle states that the actual path that a system follows will be that 
which causes (B9l to vanish. Since 77 is infinitesimal but arbitrary, the bracketed 
term in (fl9l must vanish for 55 [g] to vanish. Hence, our Euler-Lagrange equation 
is 



Pi 



' dL ' 




' dL ' 


. d a1t_ 







dL 
dq 



(50) 



For our Lagrangian in (|44l) . we have the following Euler-Lagrange equation of 
motion 



m a1t 



<9g 



(51) 



[From (K3l . we see that (BTT) is a two-endpoint equation 0.] If, for example, 
V(q) = l/2mtu 2 q 2 , (BTT) can be written as 



m^_ 2 + C*_ 2 Q - mw 2 * 







(52) 



Notice the appearance of the FFD in (BTT) . It arises because of the integration 
by parts formula ( BTT ). In order to have a strictly causal equation of motion, Riewe 
suggests considering an infinitesimal time interval [0, 2e] and then replacing all 
RFDs with LFDs. This seems unsatisfactory because fractional operators have 
memory due to their nonlocal [in time] nature. By restricting the time interval 
to an infinitesimal duration, Riewe is effectively erasing this memory. Also, it is 
questionable if this will provide an accurate approximation. For example, let our 
time period be [0, 2e] and 



/(*) 



6(t - e) 



(53) 
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Then, 

tt D*« [f] = S+Jf-e) (54) 

but, 

t F 2a [f] = y_ 2a{t -e) (55) 

Now let a = 1/2. Obviously d54l and d55l) do not agree for t < e. For t > e we 
have that ^-i(t — e) 7^ while — e) does vanish. 

If we blindly follow the above procedure for (BIT) we have the resulting equation 

dV 

mq + Cq = — — (56) 
oq 

which is missing a minus sign in front of the derivative of the potential V. We 
could of course recognize that a *F^[g] = — q for a < t < b and change the sign of 
C in (l4*4l . Then we would have the correct causal equation of motion with friction 

dV 

mq + Cq = — — (57) 
oq 

However, this requires that we treat integer ordered derivatives differently, which 
is not entirely satisfactory. 

Instead of using the Lagrangian in (14*41) . let us use 

L = -y (M [?]) (M [q]) ~ f GA" [q]) ( t D? [g]) - 7(g) (58) 
If we perturb g by 77 in d58t . we have, to first order in 77, 

8L = -y.D 4 1 [ff] t % 1 -f.D t a [g]*% a -^»7 
777 f7 1 r)V 

Then, using (BT1) . 



retarded 



fa 2 { " ^ " W } ^ (60) 



advanced 
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Now, 



5S[q] = f 

J a 



5Ldt 



(61) 



To make 8S[q] vanish, we will require that the bracketed terms in (l60b vanish 
separately. This gives us two equations of motion 



respectively, for a < t < b. 

Comparing (loH) and (l65t with ©, we see that Bateman's method is included 
in Riewe's formalism provided we use Lagrangians as in (l58l) and, require the 
advanced and retarded parts of the perturbed action to vanish separately. [These 
types of Lagrangians were not considered explicitly by Riewe in lfT2l [131 .1 Al- 
lowing both a retarded and an advanced equation of motion to arise from the 
variation of the action seems more natural than, for example, (l52t . It avoids the 
final procedure of replacing t D° [q] with a Df [q]. Also, the Lagrangian in (l58t 
is preferable to that in (l4*4l because it does not apriori assume that the LFDs are 
to be favored over the RFDs. Now we turn our attention to an alternate way of 
constructing nonconservative Lagrangians. 



In order to develop our new formalism of nonconservative Lagrangians, we will 
need some background on Volterra series J5j[l5). The Volterra series is a general- 
ization to functionals of the power series of a function. For some functional V[q], 
we define the symmetric kernels 




(retarded) 
(advanced) 



(62) 
(63) 



(64) 
(65) 



4 Volterra Series 



5 n V[q] 



(66) 



5q{r x ) ■ ■■5q{r n ) 
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The K^(-)'s are symmetric under an interchange of the r/s. So, for example, 
liQ (tljTs) = K2 (j2, Ti). Introducing the notation 

K^*q n := / Ki s \T 1 ,...,T n )q(T n )---q(T 1 )dr n ---dT 1 (67) 

we can expand the functional V[q] in the Volterra series 

V[?] = £-tf< a) *g n (68) 

71=1 n - 

[For our purposes we can assume that Kq — V[0] = 0.] It is easy to show that 

8K& * g n 



5g(t) 



/ ■■■ / r 2 ,...,r n ) (69) 

eCO • ••g(T 2 )c/r„ • ■ - dr 2 



The symmetric kernels are the natural choice to use in a Volterra series. How- 
ever, we may be given asymmetric kernels and would like to symmetrize them or 
vice versa. As motivation, consider the function 

Kq) = - 2 X ■ X A '</M, (70) 

i i<j 

where q = [qi, . . . , q n ). We can symmetrize u(q) into the form 

Kq) = I Y, K imi ( 71 ) 

where if^ = i^-, is a symmetric matrix. We will be particularly interested in 
triangular kernels given by 

K® (n, . . . , T n ) = unless n > r 2 > • • • > r n (72) 

Now, let a be a permutation of 1, . . . , n. The symmetrization of d72l is defined as 

symK® (rt,...,^) := ( 73 ) 

n! V 

= ljirW(ri,...,r B ) 
n! 
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5 Volterra Series Fractional Lagrangian Mechanics 



Let us now reconsider the nonconservative harmonic oscillator equation of motion 
in (l52l . Using the notation in (16^1 ). d52l becomes 

K^Kq 1 = (74) 

where 

Ki s \t,r) := m^_ 2 (t - r) + Ctf_ 2a (t - r) - muj 2 ^ (t - r) (75) 
Let our action be given by 

Hq] = \K { 2 s Kq 2 (76) 



Then, 



6V2[q] - (77) 



5q(t) 

Requiring (1771) to vanish gives us ( 1741 . 

Suppose now that we have a driven harmonic oscillator 

mq + muj 2 q = f(t) (78) 
We can form a new functional 

VAQ] = K^+q' + ^K^+q 2 (79) 

where 

K^(t,r) := m^_ 2 (t-r) -muj 2 ^ (t-r) (80) 
It immediately follows that, ignoring boundary conditions, 

6V2 ' [q] ~- K^(t)-mq~muj 2 q (81) 



8q(t) 



Requiring (l8*TT) to vanish and comparing with dTHl . we see that liQ, (t) = /(£) 



We can also handle higher order potentials. Let, for example, 

VAq] = q 1 + g 2 + i ^ * g 3 (82) 
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where, for some constant C, 



^f ) (n,r 2 ,r 3 ) := C^ (n - t 2 )* (t 2 - r 3 ) (83) 
Then, again ignoring boundary terms, 

M = /(t)-mg-mo; 2 g+^ 2 (84) 

We recognize <f82l> as the beginning of the Volterra series for some functional V[q] . 
To all orders of q, 

OO 1 

V[q] = £-ri4 s) *<f (85) 

n=l ""• 

[We can ignore the n = term in (l85l since this only adds an irrelevant constant 
to V[q].] For n > 2, the K^(-Ys are interpreted as the environment's reaction to 
g, which affects q's evolution. Any forcing function is included in K[ s \t). 

All of the actions considered above share two key properties: 

1. The kernels K^(-) are all localized along the line T\ = t 2 . 

2. The kernels satisfy the relation K^(-) = K^(-). 

These properties make the above actions particularly easy to analyze. However, it 
is impossible to introduce even the simple term Cq into the equations of motion 
using the \l/ a 's [see (|4*TT) 1. Using triangular, instead of symmetric, kernels results 
in a more flexible formalism. This amounts to using the $^'s in the Volterra series 
instead of the ^^'s. We will then be able to construct symmetric kernels that only 
use the $J's, not the ^q's. This requires us to be careful about the boundary terms 
in our equations. It is this situation that we now turn our attention to. 

We return again to the nonconservative harmonic oscillator. For some constant 
C, define the triangular kernels 

K+(t,r) := - [m$± 2 (t-T) +C$+ 2Q (t-r) + moo 2 ^ (t - t)] (86) 

K 2 (r,t) := --[m^z 2 (r-t)+C^Z 2a (T~t)+mco 2 ^ (r-tj\ (87) 

where 

K+{t,r) = K 2 {r,t) (88) 
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Now consider the functional 

1 rb + 



5q(t) 



V[q] := - / P ^ 2 + (r 1 ,r 2 )g(r 2 )g(r 1 )rfr 2 rfr 1 (89) 

L J a J a 

The functional derivative of (1591) is given by lfl~5l 

6V[q] {£ + / a - ir2+(ri ' T2)[g(r2)+M(r2 " t)]x 

[?( r i) + W(ri - t)]d,T 2 dT 1 - 
I [ 1 K£(T 1 ,T 2 )q(T 2 )q(T 1 )dT 2 dT 1 \ 

I 1 K 2 h (ri,r 2 )5(r 2 - t)q{T 1 )dr 2 dTi + 

J a - 

/ 1 K£(T 1 ,T 2 )q(T 2 )5(r 1 - t)dr 2 dr 1 

J a~ 

' f b ^ 2 -(t,ri)g(r 1 )tZri + ^ /* ^ 2 + (t, r 2 ) g (r 2 )dr 2 (90) 



l 

2 

l r 

2 



2 



advanced retarded 

where a < t < b. Instead of requiring the sum in d9"0l to vanish, we will require 
the advanced and retarded parts of the action's variation to vanish separately. This 
gives us two equations of motion for our system 

m$t 2 (t) + C$± 2a {t) + m^ 2 $r}(t)] * q(t) = (retarded) (91) 

m$I 2 (t) + C$Z 2a (t) + mcj 2 $o (t)] * q(t) = (advanced) (92) 

From (l9"0t. we see that g(ri) = for r x > 6 and q(r 2 ) = for r 2 < a in ( 1591 . 

Note that if our kernels only contain terms § 2n , n an integer, requiring the 
advanced and retarded parts to vanish separately is equivalent to requiring the sum 
in d9"0"l) to vanish, ignoring boundary conditions. This is because $ 2n = $ 2n an d 
both equal, up to a sign, ^> 2n . So in this case we can freely use the symmetric 
kernels ^> 2n in our action. We can also extend the above action to a driven harmonic 
oscillator and higher order potentials, as we did earlier. Again, this is due to the 
fact that = $o and also that k[ s) = Kf. 

The kernel in ( 1591 is lower triangular in the 7iT 2 -plane [i.e., K 2 (ti, t 2 ) = 
when T\ < t 2 ]. We could have equally well used the functional 

%] := I £ + f + + K+ir^rMrMn^dn (93) 
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to arrive at the equations of motion in d9"Tb and d9~2"t. Here the kernel is upper 
triangular in the r^-plane [i.e., K}(t2 : tx) = when n > r 2 ]. A derivation 
similar to that in (I^Ut shows that, if we use (1^1 for our action , then g(ri) = for 
n < a and g(r 2 ) = for r 2 > b. It follows that the symmetric action 



V[q] ■= \ {V[q} + V[q]} 



(94) 



I rb+ /-fe 4 

2 



•#2~ (ti, r 2 ) + ICj" ( r 2, n) f q(T 2 )q(Ti)dT 2 dT 



could also be used to derive (l9"Tb and d^2~l ). where g(rj) = 0, i = 1, 2, for r* < a 
and 7j > 6. The above is easier to see if we let K 2 {t) '■= K 2 (t, 0) and K2~(t) : = 
K 2 (t, 0). Then (JHU) is given by 



V[g] = - \q{t) K 2 + (t)*q(t) 



dt 



Using the integration by parts formula in (OUb gives us d93l 

1 



V[g] 



- /g(t) K 2 ~(t)*g(t) 



eft 



Adding <|95]) to (Pol) , and multiplying by 1 /2, results in $Q 



V[g] 



(95) 



(96) 



(97) 



Let us now collect some remaining observations. The usual action for the 
harmonic oscillator is given by 



S[q) 



■2 2 2 

mq — muj q 



dt 



mq + moo q dt + — qq 



(98) 



where we used an integration by parts in the second equality. The Volterra series 
action in (l95b . with C = in (l86b . gives 



V[g] 



2 
1 

2 

4 " 



m$l 9 + mw 2 $ 



mg + muj q 



mq + mu g 



* q dt 

1 
2 



g(a)5(t — a) + q(a)5(t — a) 



dt 
(99) 
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where we used (IT3t for the second equality. Hence, 



S[q] 



V[q] + - qq 



(100) 



Thus, V[q) differs from S[q) only by the boundary terms of q(t), which, by our 
above analysis, does not affect the resulting equation of motion. This also holds 
for (Pol ) and (1971 . Returning to (l95l . let us perturb q(i) by r](t). Then, 



SV[q] 



dt + 



V 



K+*q 



dt 



(101) 



advanced 



retarded 



The second term on the right of dlOll) is what we typically want in order to derive 
our equation of motion for q{t). However, the first term on the right of (1101 b is 
interesting. It shows that the advanced equation of motion for q{t) arises because 
of the perturbation of the environment's reaction due to rj{t). That is, using (l30l . 



J q [k£ *rj]dt = J n[ K 2 *q 



dt 



(102) 



So it seems that the future evolution of q(t) is affected by its past evolution because 
of the memory "stored" in the environment. 



6 Perturbing the Environment 

So far we have examined everything in a fairly standard way. We assumed that 
the environment is described by K^t) : = K^t,®) [see (ISBT ll and introduced a 
particle into this environment via q(t). Then we perturbed the particle's path and 
required that the variation in the action vanish under this perturbation. Notice 
that the particle did not change the environment's kernel given by K^t). So 
the particle itself must be so negligible that the environment's kernel does not 
substantially change under its introduction. That is, the particle is a perturbation 
to the environment. Let us explore this idea more for the nonconservative harmonic 
oscillator. This will lead to a more holistic view of mechanics which ignores the 
distinction between the environment and the system, in the case of the harmonic 
oscillator. 
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Let us assume that the environment is adequately described by the generalized 
function [see d^Ol for the second equality below] 



m$± 2 (t) + C<5>t 2a {t) + muH+it) 
m$l 2 (t) + C$± 2a (t) + m^ 2 $+(t)] * $+(t) 
^+(i)*$ + (i) 



(103) 



where the hat on $q (t) is for bookkeeping purposes only. Now let us perturb the 
environment K^t) by perturbing $q (t) by 77 (i), where fj(t) is infinitesimal in the, 
e.g., L 2 -norm compared to $q (t) = <5(t). In particular, we will not require that 
fj(t) vanish at any boundaries. Then, 



6K+(t) = K+(t)*fj(t) 



(104) 



Requiring (11041) to vanish gives us exactly (1911 ) when we identify fj(t) = q(t). 

In (11031) . we assumed that the environment reacts causally to any perturbation. 
This resulted in the retarded equation of motion in (11041) . If, instead, we considered 
the kernel K 2 (t) := K% (0, t), [see (1871) 1. then the advanced equation of motion in 
d92t would have resulted instead of (11041) . So, we see that it is not necessary, for the 
nonconservative harmonic oscillator, to construct an action in order to derive the 
equations of motion. How far this idea can be advanced to more general systems 
is an open question. 



7 Discussion 

Let us look at the retarded equation in (I9TT) a little. This is a convolution between 
the coordinate q(t) and the distribution K^t, 0) := K^t). An insightful way of 
viewing this is to think of K^t) as the environment's response "function", where, 
for an arbitrary f(t), 

y(t) = K+(t)*f(t) (105) 

Then, the actual paths that a system can follow will be those f(t) such that y(t) = 
in (11051) . That is, the path a system follows will be those such that the environment's 
response to it vanishes. This treatment of fractional derivatives as signal processors 
is well known fj"4l and can be extended to more general Volterra series than that 



19 



in (I9T1) [4, 91- Similar comments hold for (l$2"l) and its generalizations, but, in this 
case, the environment's response is anti-causal. We have not pursued this line of 
research here. However, it does open up the possibility of examining mechanics 
from a systems theoretic viewpoint [see, e.g., El for an introduction to systems 
theory] . 

We have not considered the most general action here. Instead, our attention 
was restricted to including fractional derivatives in an equation of motion derived 
by using a variational principle. In this respect we have succeeded. Further 
research is needed to see how far our formalism can be developed and, how useful 
it will be in situations other than those considered here. In particular, it would be 
interesting to extend the formalism in Sectional to more general situations. 
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